Controlled Creation of Spatial Superposition States for Single Atoms 
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We present a method for the controlled and robust generation of spatial superposition states of 
single atoms in micro-traps. Using a counter-intuitive positioning sequence for the individual poten- 
tials and appropriately chosen trapping frequencies, we show that it is possible to selectively create 
two different orthogonal superposition states, which can in turn be used for quantum information 
purposes. 



During recent years trapping and controlling small 
numbers of neutral atoms has emerged as one of the most 
active and productive frontier areas in research [H, 0, 0, Hj] . 
The interest in single atom systems is driven not only 
by the desire to perform experiments that answer long- 
standing and fundamental quantum mechanics questions 
@> 0], but also by the desire to implement concepts of 
quantum information using neutral atoms 0,0,0]. Ad- 
vances in the technology of optical lattices and micro- 
traps have recently allowed for substantial progress in 
this area 13, II , l2, EH and various concepts have been 



developed to prepare and process the states of single 
atoms in a controlled way. While techniques for control- 
ling and preparing the internal states of atoms using ap- 
propriate electromagnetic fields are well developed, only 
few concepts exist for achieving the same control over the 
spatial part of a wavefunction 14 , IB], . 



Three fundamental requirements for controlling the 
spatial part of a wavefunction are preparation, storage 
and transport. Optical and magnetic micro-potentials 
have proven to be robust tools for atom storage and tech- 
niques for moving particles between different potential 
sites using dedicated waveguides or sop histicated tun- 
neling schemes have been suggested 171 LL HH. While 
waveguides are usually highly static, the tunneling inter- 
action can be tuned by actively changing the distance 
between or the potential heights of neighboring traps. 
Both of these possibilities have recently been explored 
by Eckert et al. HI, 1(| and Greentree et al. These 
works considered three modes in three separated poten- 
tials and suggested the use of a STIRAP-like process to 
achieve a robust transfer of an atom from one trap to 
another with high fidelity. In the area of three-level op- 
tics, the STIRAP process refers to the technique whereby 
a counter-intuitive application of laser pulses leads to a 
transition of an electron between the ground states in a 
A-system 2l|, |22[ . In the atom trap scenario the energy 
levels are replaced by spatially separated trap ground 
states and the laser interaction is replaced by the co- 
herent tunneling interaction. Eckert et al. also showed 
analogues for coherent population trapping and electro- 
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FIG. 1: Left: The three atom traps are arranged in a lin- 
ear fashion. The trap on the left and the trap in the mid- 
dle are harmonic potentials and the double well trap on the 
right is made piecewise out of two harmonic traps. The trap- 
ping frequencies are arranged in such a manner that either 
the symmetric state (upper schematic) or the anti-symmetric 
state (lower schematic) of the double well is in resonance with 
the ground states of the other traps. The asymptotic ground 
states in the double well traps are indicated by a broken line. 
Right: In the counter- intuitive timing scheme the distance be- 
tween the double well and the middle trap, dMR, is reduced 
before the distance between the middle trap and the trap on 
the left, dLM- 



magnetically induced transparency [15|. They termed 
this new area three-level atom optics. 

In the following we will describe a process that ex- 
tends the existing work in three-level atom optics to the 
domain of four levels. In particular we will show how 
to create a coherent superposition of an atomic center- 
of-mass state in a controlled way and study its stability. 
While this work is based on a recent suggestion for multi- 
level optical systems 23] , its translation into the atomic 
realm does not only require considering new constraints 
on time-scales, but also allows for an easy interpretation 
of the physics involved. 

The scheme is based on two well-known, quantum me- 
chanical phenomena. The first one is the above men- 
tioned existence of dark states in multi-level systems with 
an appropriately shaped coherent coupling. In spatial 
arrangements such levels correspond to trap eigenstates 
and the coupling is facilitated by a tunneling interaction. 
No direct coupling is allowed between the initially occu- 
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pied level and the final state and this can be achieved 
by choosing a linear trap arrangement (see Fig. [1]). A 
STIRAP-like process now requires the tunneling inter- 
actions between neighboring potentials to increase and 
decrease in accordance with a counter-intuitive timing 
sequence: first the tunneling probability between the 
two empty states is increased and after a specific delay, 
At, the tunneling interaction between the occupied and 
empty state is increased. This leads to a very robust 
transfer of the particle between the two states that con- 
struct the dark state and in particular it avoids stringent 
conditions on the timing of the interactions. In a two- 
level setting such conditions are essential to avoid Rabi 
oscillations. The STIRAP process has already been ex- 
tensively analyzed in three-level optics 2lJ, |22| . 

The second phenomenon is the appearance of ground 
state splitting in double well potentials. When combin- 
ing two single traps to form a double well potential their 
respective asymptotic ground states combine to yield a 
symmetric and an antisymmetric state, with the energy 
difference between these two states depending on the dis- 
tance between the two traps. This energy difference is 
often referred to as the tunneling splitting energy and is 
directly related to the Rabi-frequency of the tunneling 
oscillations between the two wells. 

Our setup is shown schematically in Fig. [1] and con- 
sists of three traps that are arranged in a linear array 
24l | . The two leftmost traps are simple harmonic poten- 



tials, V(x) — ^muj 2 x 2 , of identical trapping frequency, 
to, and we will denote their ground states by |0)x, and 
\Q)m- The trap on the right is a double well potential 
that (for numerical simplicity) we choose to be composed 
of two simple harmonic potentials of frequency u)r, with 
asymptotic ground states 10)^ and |0)r r . The distance 
between them is fixed and given by d. The two lowest 
lying eigenstates of the double well trap are the even and 
odd combinations of these asymptotic ground states and 
are given by |0)^ = {\0)r l ± \0) Rn )/V2. Initially, a sin- 
gle atom of mass m is located in |0)l and the other traps 
are empty. After the STIRAP process has taken place, 
the atom's wavefunction will be completely transferred 
into the double well trap and possess a symmetry that 
is determined by the difference between the trapping fre- 
quencies u> — UJr. 

The time dependence of the tunneling interaction is 
realised by reducing and increasing the distance between 
the individual traps [1 51 ] - We assume the middle trap to 
be fixed and the outside traps to undergo the approach 
and reproach sequence. For the distances between the re- 
spective traps, we assume the following time dependence 

, () _jd LM iH < At 

dLM{ ' \ i [cos + 1] d* LM + d min if t > At 

(1) 

where d^M = ^lm(O) is the initial distance between the 



left hand side and the middle trap, d m i n is the min- 
imum distance between the traps in the process and 
d*LM — ^iAf(O) — dmin- We always ensure that <i m i n > 
max(a, an), where a and aR are the ground state sizes of 
the harmonic potentials for lo and wr respectively. The 
overall time for each individual approach and reproach 
sequence is T and the delay between the two sequences 
is given by At. The sequence for dMn(t) is analogous to 
eq. H]) with the constant and the cosine part reversed. 

The main advantage of the counter-intuitive timing se- 
quence is that it relaxes the stringent conditions that 
would be necessary if one tried to achieve transport by 
sequential two level tunneling. Any imprecise timing 
would lead to the occurrence of Rabi type oscillations 
and, thence, uncontrolled splitting of the wavefunction 
between all three traps. In particular, the STIRAP 
method inherently avoids the transfer of any part of the 
atom into |0)m at any time during the process. This 
is achieved by keeping the system adiabatically within a 
dark state at all times, thereby avoiding any contribution 
from the middle trap. 

To explain the dynamics when the double well poten- 
tial is involved, we will make use of a so-called double 
dark state that exists in the Hamiltonian for such a sys- 
tem. If we assume that the energy of the atom is fixed 
to the energy of the motional ground state of the left 
trap at all times, E — \TujJ, the system's Hamiltonian 
can be written in terms of asymptotic eigenstates of the 
individual harmonic traps as 
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(2) 

Here ^lLAi(t) and ^MR(t) describe the time-dependent 
tunneling frequencies between the states |0)l and \0)m 
and |0)m and \0)r , respectively. The tunneling fre- 
quency between the two states \0)r l and \0}r r is given 
by £Ir and is fixed at all times. For numerical simplic- 
ity we model the traps as harmonic oscillator potentials 
that have a fixed depth. In particular, we assume that 
this depth does not change when the traps approach each 
other, i. e. V(x) — m\n[VL(x), Vm(x), Vr(x)}. Allowing 
for the potential depth to change when the distance be- 
tween the traps changes would require a recalculation and 
readjustment of the trapping frequencies at every mo- 
ment in time, to ensure that the condition for the double 
dark state (see below) is permanently fulfilled. Recently 
it was pointed out that a situation with fixed depth traps 
can be achieved in optical traps by employing compen- 
sation lasers (l6| . 

The condition for the above Hamiltonian ([2|) to have 
an eigenstate with an eigenvalue equal to zero can be 
written as a relation between the trapping frequencies 
and the tunneling frequency within the double well trap 
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LU — LUr = ±Clj 



(3) 



Since £Ir is always positive, this condition implies that 
one dark state exists for u> > lur and one for lu < lur. 
The respective eigenstates are given by 

= cos0|O) L - sm9 [flO)^ ± \Q) Ra )/V2\ (4) 
= cos0|O) L -sm<9|0)| , (5) 



where the mixing angle 9 is defined as 
tanfl = V2- 



n 



MR 



(6) 



Note that in the case of the trap on the right hand side 
being only a single well trap , i.e. Qr = 0, this relation 
reduces to the result of 15j, i.e. all traps have to have 
identical trapping frequencies. If one fixes the distance 
between the two individual traps in the double well trap 
to be cLr = da.R, where ctR — \/Ti/ (itiujr) is the width of 
the ground state in the potential and d > 1, the tunneling 
frequency within the double well trap can be determined 
from the general relation for tunneling between two har- 
monic traps flBj 



Ofl(ef) -1 + e d [1 + d(l - crf(d))] 



,2d 2 



1)/ 2d 



(7) 



Combining this with condition ([3]), the resonance condi- 
tion for the frequency, lur, is given by 



i ± n f 



(8) 



This result can be easily interpreted. For the atom to 
move from the left to the right, the system has to sat- 
isfy energy conservation. While the asymptotic ground 
state energy levels for the double well trap, ^Ulur, are 
either larger or smaller than the ground state energies of 
the left and the middle trap, \Tilu : condition ([3]) ensures 
that either the eigenenergy of the symmetric or the an- 
tisymmetric state of the double well is equal to \Tiiu. In 
fact, for lu > lor the particle makes the transition into 
the (higher energy) antisymmetric state (^hui = E~), 
whereas for u> < lur the particle makes the transition 
into the (lower energy) symmetric state (^hio = E + )(see 
Fig.[lJ. 

To demonstrate the creation of coherent superposition 
states we have performed numerical integrations of the 
Schrodinger equation of the system, using an FFT/split- 
operator algorithm. At the start of each simulation the 
initial state is given by a Gaussian wavefunction centered 
in the trap on the left hand side and the distance between 
the two traps in the double well is given by dR = 1.5ctR. 
Typical results are shown in Fig. [2J where the time delay 
between the approach sequences is given by 10% of the 
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FIG. 2: Wavefunctions at the end of the transfer process for 
the symmetric (left column) and the anti-symmetric (right 
column) case. The panels (a) and (b) show the density and 
the panels (c) and (d) show the real (broken line) and the 
imaginary part (dotted line) of the wavefunction. The sym- 
metries are clearly distinguishable. The bottom row shows 
the change of population in the left trap (full line) and the 
two halves of the double well trap (broken line and dotted 
line). The spatial co-ordinate is scaled in units of a. 



time for each approach process, which was in turn cho- 
sen to be T = (1 .3S1jj) — 1 . The time dependence of the 
distance between the traps is chosen according to eq. ([T|) 
and the minimum distance between the traps is fixed 
at 1.2 max(a, to make sure tunneling is the only 
method for transfer. The three panels on the left corre- 
spond to the case lor > lu and the three panels on the 
right to the case lor < lu. The even and odd symmetries, 
respectively, of the final wavefunction are clearly visible. 
The two panels in the bottom row show the amount of 
|?/>| 2 in each trap. While at the beginning the wavefunc- 
tion is completely localised in the trap on the left, at 
the end of the process it is equally split between the two 
wells of the double trap. There is never any significant 
population in the middle trap (not shown). We have per- 
formed these simulations for a wide range of values for 
dR and the delay interval, At, and found this process to 
be extremely robust. 

To examine the systems fragility to noise, we simulate 
the process above using a frequency range for lur that 
is ± 10% out of resonance with the value that the sym- 
metric or the antisymmetric resonance condition ([3J de- 
mands. In Fig. [3J the upper two panels on the left (right) 
show the real and imaginary parts of the final wavefunc- 
tion for lur > lu (lur < lu). While the transferred amount 
is no longer 100% (however, it is still > 99%) for the 
non-resonant systems, one can immediately see that the 
symmetry is still a preserved property. To quantify the 
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FIG. 3: The four upper panels show the real and imaginary 
parts of the final wavefunction as a function of the trapping 
frequency in the double well trap. The trap frequency lor is 
normalised to the value required by the resonance condition 
and it is varied in an interval of ±10%. The two graphs on 
the left correspond to the case lo < lor and the two on the 
right tow> lor. The two panels in the bottom row show the 
symmetry function for the wavefunctions above as defined in 
eq. ©, where the full line represents the real parts and the 
broken line the imaginary parts. Both lines are lying on top 
of each other for the antisymmetric case. 



symmetry we define the following functions 



S R = 



E,I»W*))I 



Si = 



E^Wg)) 
EJ9#(*))I 



(9) 



which give Srj — ±1 (depending on the phase of the 
wavefunction) for a perfectly symmetric state and Srj = 
for a perfectly antisymmetric state. The lowest panel on 
the left (right) hand side shows this function for lor > lo 
{lor < lo) and confirms the optical inspection of the upper 
panels that the symmetry of the wavefunction is very 
robust against imperfections in the setup. 

The condition that the transfer has to be performed 
adiabatically is already inherent in the name STIRAP. 
The whole process therefore has to proceed slower than 
the inverse of the lowest trapping frequency to avoid ex- 
citations of higher lying levels 



T + At > 



mm la;, u r\ 



(10) 



For our setup, which includes the double trap, this con- 
dition only ensures that a full transfer will be achieved 
from the left trap into the double well trap. In order to 
achieve a steady state with the desired symmetry, a sec- 
ond condition on the timescale has to be fulfilled, which 
ensures that the process is slow with respect to the tun- 
neling frequency between the two wells in the double well 




FIG. 4: Oscillating density imbalance between the two halves 
of the double well trap, ApA = max(pz, — pr). The curve 
is for the anti-symmetric situation and plotted against the 
time interval of the process in units of j^. The inset shows 
an imbalanced wavefunction for which pL and pR describe 
the density fractions towards the left and the right from the 
center of the trap respectively. 
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If this condition is not fulfilled, the final state will no 
longer be an eigenstate but rather be given by an os- 
cillating population imbalance between the two wells. 
This, however, does not affect the symmetry of the final 
state. In Fig. 2] we show the amplitude of the oscillations 
ApA = max(pL — pr) for different values of To = T + At. 
It can be clearly seen that a steady state is established 
once condition (fTTj) is fulfilled. While the data shown 
are for the case lo > lor, the behaviour is analogous for 

LO < LOR. 

The above process shows that spatial STIRAP-like 
processes do not only allows for a robust transport or 
preparation of the wavefunctions amplitude, but also al- 
low for high fidelity when controlling the phase. It is 
therefore a good starting point for many applications in 
atom interferometry 14j. The possibility of manipulat- 



ing the parts of the wavefunction selectively with well- 
focused laser beams and closing the interferometer by 
running the whole process in reverse then opens the pos- 
sibility of creating universal quantum gates [25| . 

Schemes like this pose a challenge to current technolo- 
gies, since they require dynamical control over the po- 
sition of the micro-traps. Several different technologies 
have recently emerged that allow for such control and 
active experimental efforts are undertaken in many labo- 
ratories. Initial experiments have already shown the pos- 
sibility of dynamically controlling the distance between 
traps [H. 

In summary we have suggested a robust and straight- 
forward technique for the controlled creation of center- 
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of-mass superposition states of atoms in micro-traps. By 
adjusting the trapping frequencies the final wavefunction 
can be chosen to be symmetric or antisymmetric. This 
process relies on a STIRAP-likc method for transferring 
the atom into the new state, and the appropriate ad- 
justment of trapping frequencies in order to choose the 
symmetry. While the spatial STIRAP process has been 
shown to be very robust for transfer of amplitudes, we 
have shown that it can also be used to create a very ro- 
bust technique to control the phase of the wavefunction. 
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